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ABSTRACT 

In recent years it has become clear that understanding the absorption present in 
AGN is essential given its bearing on unification models. We present the most re- 
cent XMM-Newton observation of Mrk 6, with the goal of understanding the nature 
and origin of the complex absorption intrinsic to this source. X-ray spectral fitting 
shows that a simple warm absorption model provides an equally good statistical rep- 
resentation of the CCD data as a partial covering model. Furthermore, once the RGS 
data are included in the spectral fitting, the simple warm absorber model provides a 
very good fit to the data, without increasing the complexity of the model, in contrast 
with the partial covering model which requires the addition of either a low metalic- 
ity (<0.03 solar) thermal plasma or low temperature blackbody emission in order to 
provide a similar quality fit. The warm absorber is also a considerably more natural 
way to explain the variability observed in the X-ray absorbing column density be- 
tween the previous XMM-Newton observation and this one, requiring only a second, 
higher column density, higher ionisation, absorber to be present during the previous 
XMM-Newton observation. In comparison, the partial covering models which requires 
moving, clumpy, material relatively close to the source that result in two distinct lines 
of sight, with separate absorbing columns that each vary considerably without any 
associated change in their covering fractions, in order to explain the observed variabil- 
ity. We associate the warm absorber either with an accretion disk wind with densities 
of ^10 9 cm -3 , or with an ionised 'skin' or atmosphere of the molecular torus with 
densities of ^10 3 ^ 5 cm~ 3 . 

Key words: galaxies: active - galaxies: Seyfert - X-rays: galaxies - galaxies: Markarian 
6. 



1 INTRODUCTION 

Absorption has been widely recognised as playing a large 
role in shaping the X-ray spectra of active galactic nuclei 
(AGN), from the detailed ionised absorption features often 
observed in the soft X-ray spectrum of Seyfert 1 galaxies 
(e.g. NGC 3783; Kaspi et al. 2002, MGC-6-30-15; Fabian 
et al. 2002 ) through to the extreme cutoff of the X-ray 
continuum emission observed in Compton thick Seyfert 2's 
(e.g. NGC 1068; Bianchi et al. 2001, NGC 4945; Done et 
al. 2003). For any individual AGN, a detailed understand- 
ing of the properties of the absorption intrinsic to the source 
is essential to the understanding of the broad-band X-ray 
spectrum and, in particular, the underlying X-ray contin- 
uum emission. More generally, detailed information on the 
absorption intrinsic to a wide population of AGN gives us 
direct information on the general environment surrounding 
AGN and can test the geometry intrinsic to the standard 



unification models (Antonucci 1993, Risaliti et al. 2002). 
Observations with older X-ray spectroscopy satellites (i.e. 
Ginga, ROSAT, BeppoSAX etc) supported the implied syn- 
ergy between the optical classifications and the level of X-ray 
absorption observed in these sources (e.g. Nandra & Pounds 
1994, Smith & Done 1996). Recent work, however, has cast 
considerable doubt upon the validity of the link between 
optical type and X-ray absorption, providing the first in- 
dications that the simplistic AGN unification model may 
be significantly incorrect. In particular, the identification of 
optical type 1 sources with considerable X-ray absorption 
columns (e.g. Page et al. 2001, Akiyama et al. 2000), and the 
converse identification of optical type 2 objects with no evi- 
dence of significant X-ray absorption (e.g. Pappa et al. 2001, 
Barcons et al. 2003) , appear to be explicit contradictions of 
the implications of the standard unification model. 

It is possible to reconcile these apparent contradictions 
with the standard unification models but such explanations 
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typically invoke a specific, and often complex, absorption 
situations (e.g. Maiolino et al. 2001). Typically, the heavy, 
neutral absorption observed in type-2 AGN is comparatively 
simple; it is well represented by absorption in single medium 
composed of cold, neutral, material of varying opacity from 
N H ~10 21 -10 25+ . In contrast, some objects show consider- 
ably more complex absorption properties (both in terms of 
profile and variability). 

Studies of the physical properties of the 'intermediate' 
class of AGN, defined by Osterbrock (1981) as having weak 
composite HI line profiles (i.e. having both a weak broad and 
a narrow component) in their optical spectra, are intriguing 
because of their potential to test, and impact on, the stan- 
dard unification models. The classification of an intermedi- 
ate Seyfert galaxy (i.e. Seyfert 1.5s, 1.8s and 1.9s) is based 
on the strength of the broad component of the H/3 emission 
line relative to the broad component of the Ha emission 
line in the sense that an optical spectrum with relatively 
prominent broad H/3 is defined as a Seyfert 1.5, an optical 
spectrum with weak, but detectable, broad H/3 is defined as 
a Seyfert 1.8 and an optical spectrum with no obvious broad 
H/3 emission is classified as a Seyfert 1.9 (Osterbrock 1981, 
Cohen 1983). Typically, the ratio of the narrow components 
of the H/3 & Ha emission can be interpreted as a measure 
of the reddening of the optical line spectrum, corresponding 
to an increase in the amount of absorbing material along 
the line of sight with increasing numerical sub-class (Ward 
et al. 1987). The ratio of the broad components, by com- 
parison, is harder to interpret and, although reddening is 
one possibility, it is far from universal (see e.g. Corral et 
al. 2005). 

In the context of the standard AGN unification models, 
these line properties suggest a highly restricted range of pos- 
sible viewing angles. The viewing angle must be sufficiently 
highly inclined that the high-column density of cold material 
that forms the putative dusty molecular torus does not com- 
pletely obstruct the innermost regions of the AGN responsi- 
ble for the broad optical emission lines, yet the viewing angle 
must be small enough that there is a sufficient amount of ma- 
terial along the line of sight to significantly redden the broad 
component of the H/3 line with respect to the Ha line. Such 
a unique 'glancing' viewing angle highlights that the inter- 
mediate Seyferts are the best, and possibly the only, sources 
that allow us to study the high column density, relatively 
low ionisation, absorbing material that may form a 'skin' on 
the molecular torus or a wind from either the torus or the 
material immediately interior to it. The range of properties 
of the X-ray absorption intrinsic to intermediate Seyferts 
has proved to be considerable, both in terms of the com- 
plexity of the model required to model the absorption and 
in terms of the range of observed column densities (<10 23+ , 
Barcons et al. 2003, Schurch & Warwick 2002). Typically 
the absorption is represented either by a 'partial covering' 
model (a model in which slightly different lines of sight to 
the central continuum region sample separate absorption re- 
gions; e.g. Weaver et al. 1994, Feldmeier et al. 1999) or by 
a 'warm absorber' model (a model in which the absorbing 
material is partially ionised; e.g. Matt et al. 2004, Schurch 
& Warwick 2002). 

The Seyfert 1.5 galaxy Markarian 6 (Mrk 6), although 
well studied at longer wavelengths, was not extensively stud- 
ied in the X-ray band prior to the advent of XMM-Newton. 



Mrk 6 is an SO spiral galaxy at a distance of ~50-80 Mpc 
(Stark et al. 1992; here we use a distance of 70 Mpc) and 
was first identified as an intermediate Seyfert through opti- 
cal spectroscopy by Osterbrock & Koski (1976). A combina- 
tion of near infrared and optical data were used by Rix et 
al. (1990) to constrain the extinction, and by association, the 
column density to be A„<2 and Nh~4x10 21 cm~ 2 respec- 
tively. Mrk 6 is one of the few Seyferts to show considerable 
evidence for an ionisation cone (Meaburn et al. 1989) and 
there is some tentative evidence that the anisotropy of the 
ionising radiation is variable (Sergeev et al. 1999). MERLIN 
observations of Mrk 6 reveal collimated jets extending out of 
the plane of the galaxy to both sides of the nucleus (Kulula 
et al. 1996). These jets are remarkably similar in structure 
to those seen in Mrk 3 and are qualitatively similar to those 
seen in NGC 4151. HST observations of Mrk 6 have revealed, 
in remarkable detail, an optical emission-line jet co-spatial 
with the southern radio jet (Capetti et al. 1995). The sim- 
ilarity between Mrk 6 and NGC 4151 has been well noted 
in both the optical and X-ray regimes and Mrk 6 is often 
referred to as a "4151 analogue". NGC 4151 harbours one 
of the best studied and brightest AGN in the sky, however 
the complex behaviour and properties that are evident in its 
X-ray and optical spectra have not been observed in other 
nearby AGN. 

Despite its relative brightness, Mrk 6 has only been ob- 
served spectrally in the X-ray regime three times prior to 
the new observation presented here, once with ASCA, once 
with BeppoSAX and once with XMM-Newton. The spectral 
analysis of the 40 ks ASCA observation of Mrk 6 (Feld- 
meier et al. 1999), despite the poor signal-to-noise of the 
data, revealed a 0.5-10 keV X-ray spectrum composed of a 
canonical underlying X-ray continuum, iron Ka line emis- 
sion, and complex, heavy column density, absorption (~10 23 
cm -2 ). Feldmeier et al. interpret the heavy absorption in 
terms of a partial covering model in which the direct line 
of sight passes through a relatively dust-free torus 'atmo- 
sphere' with an several lines of sight, each of which pass 
through a significantly different column density. The X-ray 
spectrum was remarkably similar to the well-studied and 
complex Seyfert 1.5 NGC 4151, a similarity that is further 
heightened by the similar 0.5-10 keV X-ray luminosities of 
the two sources (L x ,2-io~ 5xl0 42 erg s _1 - NGC 4151; 
Schurch & Warwick 2002. c.f. L x ,2-io~l x 10 43 erg s" 1 - 
Mrk 6; Immler et al. 2003). The interpretation of Feldmeier 
et al. accounts for the complex absorption profile and re- 
solves the discrepancy between the column density measured 
by the near infrared/optical data and the X-ray data, be- 
cause the dust-free region provides photoelectric absorption 
but not optical extinction (Maccacaro et al. 1982). The re- 
cently published BeppoSAX observation yield similar details 
but with the backing of considerably better signal-to-noise 
(the LECS exposure time was ~110 ks) and, thanks to the 
large bandpass of the BeppoSAX instruments, a well defined 
and constrained continuum (Malizia et al. 2003, Immler et 
al. 2003). The BeppoSAX data also demonstrate the pres- 
ence of complex heavy X-ray absorption, which they model 
with a partial covering model similar to that of Feldmeier 
et al. (1999). Malizia et al. note, however, that the Bep- 
poSAX best-fit column densities are considerably different 
from the column densities observed with ASCA, implying 
that the absorption is considerably variable on the timescale 
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of the two observations (~2 years). Malizia et al. associate 
the measured X-ray column density with absorption origi- 
nating in the broad line region, however it is unclear from 
the data whether the absorbing medium requires a partial- 
covering model, a warm absorption model, or indeed a com- 
bination of the two. No significant variability was detected 
in either the ASCA or BeppoSAX observations of Mrk 6. 
BeppoSAX appears to have observed Mrk 6 in a brighter 
state than ASCA, resulting in a somewhat different over- 
all spectral shape. The first XMM-Newton observation of 
Mrk 6, in conjunction with the BeppoSAX data, showed a 
similar situation, highlighting the similarities between Mrk 
6 and NGC 4151 and revealing long timescale absorption 
variability (Immler et al. 2003). Despite having ~7 times the 
counts of the previous ASCA data, the XMM-Newton data 
were still insufficient to clearly distinguish between a partial- 
covering model, a warm absorption model, or a combination 
of the two for the intrinsic absorption. The ability to distin- 
guish between these scenarios was hampered in particular 
by the poor quality of the soft X-ray spectrum both in the 
EPIC and RGS instruments. 

Here we present the most recent XMM-Newton observa- 
tion of Mrk 6. The observation was performed with the aim 
of understanding both the nature of the soft X-ray emission 
through the RGS data and the X-ray absorption through 
the EPIC instruments. This paper is organised as follows. 
Section 2 details the XMM-Newton observation. Section 3 
presents the spectra and discusses the spectral modelling of 
both the soft X-ray emission and the absorption. Section 4 
compares the details of the current observation to the pre- 
vious data and the implications that arise as a result. 



2 THE XMM-Newton OBSERVATIONS 

Mrk 6 was observed with XMM-Newton on April 26 J \ 2003 
(orbit 619) for a total of 55 ks. The EPIC and RGS in- 
struments were active for the full duration of the observa- 
tion however considerable prolonged background flaring was 
present during the first 10 ks of the observation and as a 
result the EPIC cameras were operated in closed mode dur- 
ing this period. For the remaining 45 ks of the observation 
the EPIC MOS and PN operated in Full Window Mode 
with the Medium filter in place. The data were processed 
with the standard EPIC and RGS processing chains incorpo- 
rated in the XMM-Newton SAS v6.1.0. The EPIC data were 
processed using the latest CCD gain values and only X-ray 
events corresponding to patterns 0-12 in the MOS cameras 
and 0-4 in the PN camera were accepted. An examination 
of the full field lightcurve (excluding the emission associated 
with Mrk 6 and other identified point sources) revealed sev- 
eral further periods of brief, but considerable, background 
flaring. Periods during which the background contribution to 
the count rate within a 2' radius region centred on Mrk 6 ex- 
ceeded 10% of the source countrate in this region have been 
excluded from the subsequent analysis, resulting in total ef- 
fective exposure times of 34 ks and 28 ks for the MOS and 
PN instruments respectively. A brief analysis of the impact 
of pile-up on the observation showed that the effect was neg- 
ligible due to the moderate brightness of the source. EPIC 
lightcurves and spectra of Mrk 6 were extracted from a 1' ra- 
dius circle centred on the source. Background lightcurves 
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Figure 1. The EPIC MOS (Black - MOS 1, Red - MOS 2) and 
PN (Green) spectra of Mrk 6. The lower panel shows the ratio of 
the data to a basic absorbed power-law fit. 



and spectra were extracted from source free regions of the 
central CCD's of the MOS detectors. The background sub- 
tracted EPIC MOS lightcurve is consistent with being con- 
stant for the duration of the observation (x 2 =445 for 410 
d.o.f) with an average combined MOS count rate of ~1.7 
c/s. 

The RGS instruments were active for the full 55 ks du- 
ration of the observation but are subject to the same flaring 
periods as those observed in the EPIC cameras which, again, 
have been excluded from the the subsequent analysis. Dis- 
persed RGS source and background spectra were extracted 
from the quiescent periods with the automated RGS extrac- 
tion tasks in the SAS. The 2 nd order RGS spectra contain 
no useful data, primarily because of the relatively faint na- 
ture of the source, and as a result only the 1 st order spectra 
are presented here. The resulting RGS spectra have a total 
effective exposure time of 34 ks for both the RGS 1 and 2 
instruments. 

Spectral fitting was carried out using the XSPEC vll.3 
software package with instrument response matrices gener- 
ated by the SAS software for each of the spectra indepen- 
dently. The source spectra were grouped to a minimum of 
20 counts per bin to allow the use of \ 2 minimisation tech- 
niques. 



3 MODELLING THE XMM-Newton SPECTRA 
3.1 The EPIC Spectra 

We begin by fitting the EPIC MOS and PN data between 
1-5 keV and 7.5-10 keV with a basic absorbed power-law 
model. This model, while obviously unrealistic, allows us to 
characterise the basic properties of the X-ray spectrum, par- 
ticularly for comparison with other AGN. Figure 1 shows the 
EPIC MOS and PN spectra of Mrk 6, along with the ratio of 
the data to a basic absorbed power-law fit. The fit is charac- 
terised by a very flat (F~1.2) continuum and absorption an 
order of magnitude above the level of Galactic absorption in 
the direction of Mrk 6 (7xl0 21 cm" 2 c.f. 6.4 xlO 20 cm" 2 - 
Stark et al. 1992). The positive residuals around 6.4 keV are 
indicative of the presence of a neutral iron Kq line and the 
large positive residuals below 0.8 keV demonstrate the pres- 
ence of a strong 'soft X-ray excess', a feature that is common 
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to many AGN (e.g. Piro, Matt & Ricci - 1997, Gierliriski & 
Done - 2004). The extremely flat continuum in this simple 
fit is suggestive of a strong Compton reflection component, 
coupled with a steeper underlying continuum. Indeed, anal- 
ysis of the high energy (>10 keV) X-ray spectrum of Mrk 
6 from a previous BeppoSAX observation confirms this, de- 
tecting a steeper underlying continuum (1.6<F<2.0) and a 
Compton reflection component with a reflection fraction of 
~1.0 (Malizia et al. 2003; Immler et al. 2003). 

Both the MOS and PN spectra show similar large-scale 
features with one exception highlighted by the fit to the 
absorbed power-law model. Between 0.5-0.6 keV, the PN 
shows a considerable difference in spectral shape to both 
MOS spectra which both show a dip in their spectra. An 
investigation of the RGS 1 spectra supports the MOS spec- 
tral shape in this region (the RGS 2 spectrum suffers from a 
defective CCD across this energy range). We initially limit 
our modelling to the 0.8-12 keV region for both the MOS 
and PN cameras to exclude the soft X-ray excess and sim- 
plify the analysis of the continuum, given the discrepancy 
between the PN and the other instruments we retain the 0.8 
keV lower limit for the PN data in the broadband spectral 
fits detailed in Section 3.3. 

We define a base model, based on the phenomenological 
properties observed in the over-simplistic absorbed power- 
law fit, with the following components: 

• A power-law continuum. The continuum slope is ini- 
tially fixed to be r=1.81, as determined from the analysis 
of the high energy BeppoSAX data by Immler et al. (2003), 
while the normalisation, Ai, remains a free parameter in the 
fits. 

• A neutral Compton reflection component (modelled by 
PEXRAV in XSPEC, Magdziarz & Zdziarski 1995). The re- 
flection scaling factor, R, is fixed to be 1.21, the best-fit 
value found from the analysis of the BeppoSAX data by 
Immler et al. (2003) as are the other remaining param- 
eters (i.e. incident continuum slope, r=1.81; Normalisa- 
tion, A2=4.68xl0 -3 photon cm' 2 s _1 keV -1 ; Inclination 
cos z=0.94). In addition the metal abundance in the reflec- 
tor was fixed at the solar value. 

• Complex absorption. Two distinct models were used for 
the nature and form of the complex absorption. The first 
model is the partial covering model favoured in previous 
analyses of the X-ray spectrum of Mrk 6. The second model 
is a layered model that incorporates a heavy column density 
(e.g. ~10 23 cm -2 ) partially ionised absorber, a model that 
has been extremely successful in modelling the detailed X- 
ray spectrum of NGC 4151. The details of the absorption 
modelling are the primary focus of the following discussion 
and are given in detail there. 

• An intrinsically narrow iron Ka emission line of inten- 
sity Ikc , at an energy E Ka ■ 

• Galactic absorption (6.4xl0 20 cm -2 ). 



3.1.1 The partial covering model 

As discussed in Section 1, previous approaches to modelling 
the absorption intrinsic to Mrk 6 have been largely based on 
a partial covering model (e.g. Feldmeier et al. 1999; Immler 
et al. 2003). In partial covering models, a fraction, f CO v, of 
the hard power-law continuum is absorbed by a cold gas with 
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Figure 2. The EPIC spectra of Mrk 6 fit with the base model 
incorporating partial covering absorption. 



column density Nh,i, whilst the remaining fraction of the 
continuum, 1 — fcov, is absorbed by a considerably different 
column density Nh,2- This absorption profile was incorpo- 
rated into the base model defined previously and fit to the 
three EPIC datasets simultaneously, resulting in a reason- 
able fit (x 2 =2337 for 2263 d.o.f; Figure 2; Table 1 - Model 
PCI). Despite the fit being acceptable, there is a clear de- 
viation of the model from the data above energies of ~8 
keV. The deviation is particularly evident in the PN data, 
but is also observed to a lesser degree in the MOS datasets. 
The high energy deviation is suggestive of a mismodelling 
of the Compton reflection continuum which, given the tight 
constraints employed in the base model, is not surprising. 
Freeing up the relative reflection component, R, results in a 
considerable improvement to the fit (Ax 2 =176 for 1 d.o.f; 
Table 1 - Model PC2) with a reflection component ~3.3 
times stronger than the previous measurement. Despite the 
considerable improvement of the fit the PN data >9 keV 
still lie systematically above the model. It is not clear that 
the PN calibration >10 keV is well understood, however 
under the assumption that it is, the under-prediction of the 
model thus suggests that the continuum is incorrectly mod- 
elled by the tightly constrained power-law model. A careful 
investigation of x 2 - s P ace when the photon index of the un- 
derlying continuum, F, is allowed to be a free component in 
the fit reveals a considerably better fit to the data (Ax 2 =44 
for 1 d.o.f; Table 1 - Model PC3) with a resulting best-fit 
power- law photon index of ~1.44 and a considerably weaker 
contribution from Compton reflection (R~0.6). The surpris- 
ingly flat slope appears to be well constrained, however the 
corresponding change in the strength of the inherently flat 
Compton reflection indicates a strong degeneracy between 
the two parameters shown clearly by Figure 3. 

We reiterate that the discrepancy between the model 
and the data was only present in the PN data for the pre- 
vious fits and it is possible that a calibration difference 
is responsible for the over-prediction. Allowing the relative 
reflection component and the continuum photon index to 
take different values for the fit to the PN data than for 
the MOS data to test this possibility resulted in two con- 
siderably different best-fit values for the continuum photon 
index and Compton reflection. The best-fit MOS values are 
r=1.40+°;j2, R<0.27, in contrast to the best-fit PN values of 
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Table 1. The best-fitting EPIC spectral parameters 

Parameter PCI PC2 



PC3 



WA1 



WA2 



fco V , 1 



0.57 
0.43 



+0.03 
-0.03 
+0.05 



0.44 

0.37 



+0.05 
-0.05 
+0.05 



0.28 
0.36 



+0.04 
0.03 
+0.03 
0.04 



0.47 1 



0.028° 



log£ 

fcov ,2 



3.95 
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+0.35 



2.22 
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+0.26 



1.39 



+0.0G 
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1. 
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1.21 c 
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'-0.4 



1.81 
4.74 
4.00 
6.44 
1.0 



+0.42 
0.44 
+0.32 
0.32 
+0.03 
0.03 
+0.4 
0.3 



1.44 
3.10 
0.60 
6.44 
1.34 



+0.09 
0.07 
+0.20 
0.21 
+0.14 
0.14 
+0.03 
0.02 
+ 0.36 
0.34 



1.81 c 



6.18 



+0.07 



6.45 



0.07 
1.21 c 
+0.02 

0.03 



1.62 



+0.33 
0.34 



1.81 

4.67 
4.3 
6.44 
0.93 



+0.10 
0.08 
+0.2 
0.3 
+0.04 
0.03 
+0.37 
0.34 
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Figure 3. Confidence contour plot of T and R for Model PC3. 
The black, red & green lines are confidence contours corresponding 
to 1, 2 and 3<r, respectively (Ax 2 =2.71, 4.0 & 9.0). The diagonal 
nature of the contours highlights the strong degeneracy between 
the parameters. 

r=1.57lo ol, R=1.0l° |, again highlighting the degeneracy 
between photon index and Compton reflection. 

We stress that the fits with T~1.4 and T>1.6 both pro- 
vide a statistically acceptable fit to the data and, given that 
such a flat continuum slope is incompatible with the previous 
wide bandpass BeppoSAX data (which in several analyses 
has measured a best-fit photon index of F>1.6 using similar 
models to that presented here), we argue that Model PC2, 
despite a high contribution from Compton reflection and 
residuals in the PN fit >9 keV, is a more plausible physical 
representation of the data. 

3.1.2 The layered model 

The second approach to modelling the absorption intrinsic 
to Mrk 6 utilises a model in which a heavy column density 
partially ionised absorber is layered with a lower column 
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Figure 4. The Data/Model ratio residuals resulting from fitting 
the EPIC spectra of Mrk 6 with the base model incorporating 
partially ionised absorption. 

density neutral absorbing column. This model has been ex- 
tremely successful in detailing the X-ray spectrum of NGC 
4151 (Schurch & Warwick 2002, Schurch et al. 2003). In this 
model the absorption affects all the power-law continuum 
flux, instead introducing the required complexity in the ab- 
sorption profile through the details of the ionised absorption. 
Here, we use the same warm absorption model, calculated 
by XSTAR, that was used in the Schurch et al. (2003) anal- 
ysis of NGC 4151. This absorption profile was incorporated 
into the base model defined previously and fit to the three 
EPIC datasets simultaneously, resulting in an acceptable, if 
marginally worse, fit to that of the initial tightly constrained 
partial covering model (x 2 =2365 for 2264 d.o.f; Fi gure 4; Ta- 
ble 1 - Model WA1). Clearly the fit suffers from the same 
deviations from the data at high energies that were observed 
with the partial covering model fit, along with somewhat 
more pronounced residuals <1.5 keV. The change in the 
model of the absorption profile does not strongly impact on 
the model fit >5 keV, as highlighted by the reaction of this 
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fit when the relative reflection parameter is allowed to be a 
free parameter (Table 1 - Model WA2); the fit improves con- 
siderably (Ax 2 =89 for 1 d.o.f) and the reflection component 
is ~3.6 times stronger, a similar effect to that previously ob- 
served with the partial covering model. The increase in the 
strength of the reflection component is coupled with con- 
siderable decreases in both the ionisation state of the warm 
absorber and its column density. Similar effects to those de- 
tailed in the previous section are also observed when the 
photon index, F, is available as a free parameter in the fit 
and when the parameter values for the PN data are decou- 
pled from the values for the MOS data. 

3.2 The RGS spectra 

The only appreciable difference between the best-fits of the 
two absorption models to the EPIC data is the presence 
of more significant residuals in the 0.8-1.5 keV band in the 
warm absorber best-fit. Given the apparent similarity be- 
tween Mrk 6 and NGC 4151 it is not unreasonable to suspect 
that the soft X-ray emission in Mrk 6 is composed almost 
entirely of discrete line emission similar to that in NGC 4151 
(Schurch et al. 2004, Ogle et al. 2000). In such a scenario, 
the presence of increased residuals around 0.8-1.5 keV can 
be directly attributed to the presence of this line emission 
(specifically, Magnesium and Neon emission lines). To inves- 
tigate this possibility in detail it is prudent to investigate the 
RGS data, with a view to identifying specific emission line 
features. 

The RGS data are shown in Figure 5 along with the 
ratio of the data to a basic absorbed power-law fit, in a sim- 
ilar fashion to Figure 1. The RGS spectra clearly demon- 
strate that the soft X-ray emission in Mrk 6 is dominated 
by some form of continuum emission, very different in na- 
ture to the soft X-ray line emission in NGC 4151. Despite 
this difference, the emission does show features consistent 
with possible line emission or edge absorption, particularly 
in the 0.4-0.7 keV and 1.3-1.9 keV regions, modelling of the 
RGS data independently of the EPIC data reveals that the 
soft X-ray spectrum is not well fit with simple models in- 
cluding a thermal plasma (x 2 =280 for 95 d.o.f), emission 
from a blackbody (x 2 =H5 for 95 d.o.f) or a simple absorbed 
power-law (x 2 =179 for 96 d.o.f). Following these simple fits 
we fit the RGS spectra with the most reasonable fits from 
the EPIC analysis (Models PC2 and WA2). Most of the free 
parameters in the RGS fits were fixed to the values from the 
best-fits to the EPIC data, either because the components of 
the model they refer to are outside the energy range of the 
RGS (e.g. the iron Ka line) or the components are broad- 
band continuum components that are considerably better 
constrained by the EPIC data. In the best-fit partial cover- 
ing model (PC2) , the power-law normalisations and the low- 
est column density absorber were allowed to be free parame- 
ters in the fit. The model provides a barely reasonable fit to 
the RGS data (x 2 = 130 for 95 d.o.f); in addition, the column 
density of the absorption drops to ~2x 10 21 cm -2 , well below 
the best-fit value from the EPIC data, and the model leaves 
considerable residuals around the relatively complex 0.4-0.7 
keV and 1.3-1.9 keV regions. Adding a thermal plasma com- 
ponent (MEKAL) does not significantly improve the fit or 
allow the low column density to take on a more reason- 
able value. Allowing the chemical abundance of the thermal 
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Figure 5. The RGS (Black - RGS 1, Red - RGS 2) spectra of 
Mrk 6. The lower panel shows the ratio of the data to the basic 
absorbed power-law fit in the upper panel. The gaps in the data 
correspond to the energy ranges covered by the failed CCDs in 
the RGS detectors (RGS1: ~0.8-1.25 keV, RGS2: ~0.5-0.66 keV) 



Table 2. The best-fitting RGS spectral parameters 



Parameter 


PC2+Blackbody 


WA2 


lvr a 
1N H,warm 

log? 

A p; c 
kT d 

Aw, 6 


i-04«; 22 

3.6±° :i 

15+ 004 
u - lo -0.03 

92+ 16 


n 99 +0.05 
U - ZZ -0.11 

2-311:1 

1 93+0.17 
4 63+1.06 


x 2 


76 


83 


d.o.f 


93 


93 


a XlO 22 an" 2 






^ Upper Limit 






10 ~ 3 photon 


keV -1 cm~ 2 s _1 




d keV- 1 






e 10~ 5 photon 


keV -1 cm~ 2 s _1 





plasma to be a free parameter in the fit does improve the fit 
but results in a chemical abundance <3% of solar abundance 
values, essentially removing all of the line emission from the 
model. Adding a low temperature blackbody component in 
place of the thermal plasma provides a considerable increase 
in the quality of the fit (Table 2) and may represent the tail 
of the termal X-ray emission from a the accretion disk. We 
note that the properties of the black body component in the 
warm absorber model are poorly constrained due to an obvi- 
ous degeneracy between the normalisation of the power-law 
continuum and the temperature of the black body emission. 

By contrast the best-fit EPIC warm absorption model 
(WA2) provides an excellent representation of the RGS data 
(x 2 =83 for 93 d.o.f) without the addition of new spectral 
components. In this fit both the warm and cold column den- 
sities, along with the power-law normalisation, remain free 
parameters. The fit differs from the best-fit EPIC values, 
suggesting a slightly more ionised medium, a smaller col- 
umn density of warm material and considerably more cold 
material long the line of sight (Table 2). The considerably 
larger cold column density is required to prevent the warm 
absorber model over-predicting the soft X-ray emission. 
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Figure 6. The EPIC and RGS spectra of Mrk 6 fit with the best- 
fit partial covering absorption model including emission from a 
blackbody. 



3.3 Joint, broadband, spectral fitting 

The RGS spectra for Mrk 6 highlight the importance of fit- 
ting the full energy range of the spectra when using warm ab- 
sorption models unless there is unambiguous evidence that 
the leaked flux from such a component is not present in the 
soft X-ray emission (as is the case in NGC 4151). Taking 
into account the details of the best-fits to the RGS data, 
we performed joint spectral fits on the EPIC and RGS data 
over the full 0.3-12 keV energy range for the MOS detec- 
tors and a more limited energy range for the PN detector 
(the PN data up to 0.8 keV remained excluded for the rea- 
sons discussed in Section 3.1). The initial models fit to the 
joint EPIC and RGS spectra have the same free parame- 
ters listed earlier for the WA2 and PC2 models however, in 
keeping with the RGS fits, a simple blackbody model was 
added to the PC2 model. Figures 6 & 7 show the best-fits 
of the partial covering absorption model and the warm ab- 
sorber model to the EPIC and RGS data, respectively. Both 
models are a good statistical representation of the data with 
reasonable physical parameters (x 2 =2327 and 2360 for 2419 
and 2422 d.o.f, respectively). Details of the fits are given in 
Table 3. The warm absorber fit is somewhat improved by the 
addition of a scattered power-law component (Ax 2 =37 for 
1 d.o.f) however the statistically good nature of the fit prior 
to the addition of the scattered power-law prevents us from 
concluding that this component is required in the modelling. 
When included, the scattered power-law contains ~3.5% of 
the flux of the underlying continuum. 
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Figure 7. The Data/Model ratio residuals resulting from fitting 
the EPIC and RGS spectra of Mrk 6 with the best-fit warm ab- 
sorption model. 



Table 3. The best- fit parameters for the joint EPIC and RGS 
spectra 



Parameter 


PC2+bb 


WA2 




0.891°;" 


27+ - 1 


fcov,l 








□ 8+1.3 


2 89+°- 14 


log? 




2 069+ 025 
z.uo»_ 021 


fcov^l 


U - OJ --0.08 




r p ; 


1.81 c 


1.81 c 


A pl d 


4 48+ - 78 


4 6+ 01 


R 


3 47 +0.34 
0.38 


4 44+0-26 
0.26 


Ejf a 


fi 44+0-03 
D ' 44 -0.03 


O-+4-0.03 


l Ka e 


0.98i°j3 


91+0- 35 

"•"—0.34 


kT d 


o.i7± ; l 




Am,' 


0-6812:1? 




x 2 


2327 


2360 


d.o.f 


2419 


2422 



a X10 22 cm~ 2 

fa Upper Limit. 

c Fixed parameter. 

d 10~ 3 photon keV -1 cm -2 s -1 

e 10 J photon cm~ 2 s — 1 

f 10~ 5 photon keV -1 cm~ 2 s _1 



4 DISCUSSION 

The analysis presented in Section 3 demonstrates that the 
absorption profile of Mrk 6 is extremely complex and that 
distinguishing between the plausible absorption profiles is 
a tricky endeavour. Despite the partial coverer providing a 
marginally better broadband fit to the XMM-Neuiton data, 
both the partial covering model and the warm absorption 
model are good representations of the data with reasonable 
physical parameters and there is little to distinguish between 
the two. 



4.1 A comparison with previous XMM-Newton data 
(1997) 

Any model of the absorption in Mrk 6 needs to be both a 
good statistical representation of the data in any specific 
observation and needs to be able to explain any observed 
absorption variability in a physically plausible manner. We 
stress the latter clause because a comparison of the results 
from the spectral fitting from the observation presented here 
to the previous XMM-Newton observation clearly indicates 
that considerable Nh changes are evident in Mrk 6, regard- 
less of which absorption model you use to represent the data. 
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4-1.1 Partial covering 

Analysis of the previous X-ray spectra from Mrk 6 have 
favoured a 'duel partial covering' model, the details of which 
are considerably more complex than the partial covering 
model presented here (Feldmeier et al. 1999, Immler et 
al. 2003 & Malizia et al. 2003). The model used in Section 3.1 
defines two distinct lines of sight to the source, both with 
different intrinsic amounts of absorption. Such a model is 
only plausible in the situation where the absorbing medium 
is sufficiently close to the continuum source that the con- 
tinuum source does not appear point-like with respect to 
the absorbing material and that the region over which the 
column density changes significantly is small (i.e. that the 
edges of absorbing clouds are relatively sharp). Taking the 
assumption that the X-ray continuum emission originates in 
the inner regions of the accretion disk, this would imply that 
the absorption originates in the broad-line region clouds, or 
closer material (for a 10 M^black hole, an emission region 
of 5R S corresponds to ~3'at a distance of 0.1 pc, compared 
to 20"at a distance of lpc). 

In fact, the 'duel partial covering' models favoured in 
previous modelling of Mrk 6 define a total of four unique 
lines of sight, three with different, distinct, intrinsic absorp- 
tion and one intrinsically unabsorbed line of sight. Despite 
the good statistical fit found with a considerably less com- 
plex model, involving only two distinct lines of sight, we fit- 
ted the XMM-Newton spectra with the more complex model 
for the purposes of direct comparisons with the previous 
studies. This model (Table 4 - XMM-Newton 2003) results 
in fit of similar quality to the initial partial covering model 
(X 2 =2330 for 2363 d.o.f). Examining the details of the four 
different lines of sight in this model, it becomes clear that 
two of the absorbed lines of sight are strongly dominant, ac- 
counting for a total of 50% and 47% of the continuum flux 
and with column densities of ~4.5xl0 22 and ~0.63xl0 22 
cm~ 2 , respectively. These values are extremely similar to 
the values for the two absorbed lines of sight used in the ini- 
tial partial covering model (c./. Table 1 - Model PCI) con- 
firming that despite greater complexity the 'double partial 
covering' model closely resembles the simpler fit. Similarly, 
adding an unabsorbed continuum component to the simple 
initial partial covering model does not improve the spec- 
tral fit considerably (Ax 2 = l for 1 d.o.f) and results only in 
an upper limit of 3.5xl0~ 5 photon cm~ 2 s _1 keV -1 on the 
flux in the unabsorbed component (<0.3% of the flux in the 
absorbed components), again confirming that there is little 
additional complexity in the X-ray spectrum >0.8 keV. 

A direct comparison between the best-fit column den- 
sity values for the 2001 and 2003 observations (Table 4) 
reveals a factor of ~2-3.5 drop in the column density dur- 
ing the intervening period. Remarkably, despite such a large 
column density variation, the partial covering fits still re- 
quire several lines of sight to the source and, furthermore, 
the best-fit lines of sight are very similar to those present 
in the earlier observation. This argues strongly against the 
partial covering model being the correct interpretation of 
the complex absorption since it is difficult to conceive of a 
situation in which absorption in clouds close to the source 
results in large column density variations without consider- 
able differences in the covering fractions of the clouds. 



Table 4. The 'dual partial covering' model 

Parameter XMM-Newton 2001 XMM-Newton 2003 



N H ,i° 

fcov,l 
fcov,2 
fcov,3 
fcov 



J- pi 

R 



X 

d.o.f 



10 ^2 +(:i ' 72 

u - ol -0.02 
S 05+ ' 69 
-+0.01 
-0.01 

2.27+ 03 



4.43 
0.50 



0.05Z 



0.40 



0.08 
+0.01 



+0.17 
0.33 
+0.01 
-0.01 

3.8+"-" 

— 0.35 

02+ 01 
u - uz -o.oi 

,+0.01 
0.01 
+ 0.01 



0.63 
0.46 



+0.01 



0.04 



1.81 c 

c oq+0.05 

1 21+ ' 14 
^^-o.ie 

6 47+0-03 
—0.03 
d+0.39 
0.53 



0.02 



+0.01 



1.68: 



5.98 



6.45 
1.7 



1.81 c 

+0.04 
0.03 
1.21 c 
+0.02 
0.03 
+0.3 
0.4 



813 
885 



2330 
2263 



XlO^ cm"' 
Unabsorbed line of sight 
Fixed parameter. 
10~ 3 photon keV -1 cm~ 2 s~ 
10~ 5 photon cm~ 2 s _1 



4-1.2 Warm absorption 

Despite not being the favoured model for previous obser- 
vations of Mrk 6, the warm absorption model does provide 
a good fit to the XMM-Newton data. Comparing the de- 
tails of the current warm absorber model parameters to 
those given in Immler et al. (2003) for the 2001 observa- 
tion reveals considerable changes in both ionisation and col- 
umn density of the absorption (log£/ i=2.42, log£/ 03 =2.07; 
N H /oi = l-33xl0 23 cm~ 2 , N ff/03 =0.29x 10 23 cm -2 ). Assum- 
ing that the absorption originates in a single 'cloud' (i.e. 
a region of material with a single ionisation state and den- 
sity), these differences imply very different properties for the 
cloud seen in the 2001 observation than for the cloud present 
during the 2003 observation. In particular, the considerable 
change in ionisation state implies either a large change in 
the radial distance, or a large change in the cloud density. 

An alternative interpretation can be made in terms of 
two separate ionised regions, a uniform low column density, 
low ionisation region and an irregular higher column den- 
sity higher ionisation region. There is considerable evidence 
for absorption with multiple ionisation states in Seyferts, 
although as yet no multiple warm absorbers have been iden- 
tified with such a heavy column density. In this model, the 
lower column density, lower ionisation state region remains 
constant between both observations, while the second re- 
gion (comprising a higher column density, higher ionisation 
state) is only present along the line of sight during the first 
XMM-Newton observation and has a higher ionisation state 
either as a result of a closer proximity to the ionising con- 
tinuum emission or as a result of a lower density. We note 
that this is not the same as the partial covering model since 
it requires only one line of sight. 

Fitting the EPIC PN and MOS spectra from the 1997 
observations with this multi-zone warm absorber model 
(Figure 8 & Table 5) results in an excellent fit to these data 
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Figure 8. The 2001 EPIC spectra of Mrk 6 fit with the best-fit 
multi-zone warm absorber model. 



Table 5. The best- fit parameters for the XMM-Newton 2001 spec- 
tra. 

Parameter Multi-zone warm absorber 



N fl ,i" 
log£i 

^H.cold" 

A pl c 
R 

A ^ 

x 2 

d.o.f 



11.61& 

o ci O+0.031 
z.oj.^_q 052 



2.89 
2.069 



+0.25 
0.35 
1.81 6 
+0.16 
0.34 
q+0.52 



1.25 
1.1 

4.39 

1 9"?~ t 
±,1,o -0.59 
i q+0.07 



6.44 



1.47 



+0.38 



1866 
2425 



a XlO 22 cm" 2 

Fixed parameter. 
c 10 -3 photon keV -1 cm -2 s _1 
d 10~ 4 photon keV -1 cm -2 s _1 
e 10 ~^ photon cm -2 s — 1 



with reasonable physical parameters and, with the heavier 
warm absorbing material no longer present along the line 
of sight, also results in a good fit to the current XMM- 
Newton data. 

We note that the model still includes a significant 
amount of cold absorption as a free parameter in the fit 
and a scattered power-law (similar to that detailed in Sec- 
tion 3.3). There is a considerable decrease in the cold column 
density between 2001 and 2003 implying that it is not asso- 
ciated with the large scale galactic environment of Mrk 6. 
The drop in the cold column density results in the scattered 
power-law becoming an apparently more important compo- 
nent in the model fit to the 2003 observations. Despite this, 
the scattered power-law only contains ~2.5% of the flux of 
the underlying continuum in the fits to the 2001 data, an 
extremely similar value to that found from the 2003 XMM- 
Newton observations. 



4.2 The origin of the complex absorption in Mrk 6 

Despite the progress made on the characterisation of the 
absorption in Mrk 6 it remains difficult to make any con- 
crete statement on the physical origin of the absorption. 
The two models discussed here imply very different origins. 
The partial covering model suggests that the origin of the 
absorption arises from material in the broad line region, or 
closer, and that the absorbing column changes considerably 
(presumably through clouds moving across the line of sight) 
without significantly changing the partial covering fractions. 
Given the complexity of the partial covering model as it is, 
it is not possible to draw any clear conclusion on whether or 
not the partial covering material is genuinely cold from the 
data. However, under the assumption the partial covering 
material is not significantly ionised then the material must 
reside a considerable distance from the ionising source (the 
absorption profile from material with an ionisation param- 
eter of <1 is not considerably different from that of strictly 
cold material, and given an ionising luminosity of 10 43 erg 
s" 1 and a density of 10 12 cm -3 results in an a cloud dis- 
tance >100R S ). This does not rule out absorption in a wind 
from the accretion disk, particularly given the lack of tight 
constraints on the possible ionisation states of the partial 
covering media. 

Similarly, possible origins for the warm absorber inter- 
pretation include an ionised skin on the torus, warm gas in 
the vicinity of the broad line clouds or a disk wind, however 
in this instance the ionisation state can give us a handle 
on how plausible an origin in each of these regions is. For 
a given distance, ionising luminosity and measured ionisa- 
tion parameter we can calculate the density of the absorbing 
material. Examination of these density values in compari- 
son to the densities expected for the material in these re- 
gions can help us locate the absorbing material (we assume 
a 10 s Mq black hole, an ionising luminosity of 10 43 erg s _1 
and distances of 60R S , 10 light-days and 0.5pc for the disk 
wind, BLR and the inner edge of the torus respectively). Lo- 
cating the absorption in a torus skin results in densities of 



^10 d 



cm which, at the upper limit, is consistent with 



the values calculated by Risaliti, Elvis and Nicastro (2002) 
for clouds in the torus with column densities of the order 
of 10 23 cm -2 . Locating the ionised material in the BLR 
(at ~10 light-days; Peterson 1993) results in densities of 
~10 6 ^ 7 cm -3 , considerably lower than the densities inferred 
for several AGN from ultraviolet observations (e.g. Gaskell 
& Sparke 1986). Placing the absorption in a disk wind close 
to the ionising source results in densities of ~10 9 cm -3 which 
are in excellent agreement with the disk wind simulations 
by Proga & Kallman (2004). Given the somewhat low val- 
ues implied for a location at distances appropriate with the 
torus and BLR we favour a location in a disk wind extend- 
ing on scales of tens to hundreds of Schwartzchild radii. In 
such a scenario the 'double warm absorber' seen in the 2001 
XMM-Newton observation is easily explained as a line of 
sight which passes through two regions of the wind, one 
located somewhat closer to the ionising source (and hence 
more ionised) than the other. The apparent variability of 
the column density is then easily explained as an effect of 
the clumpy, outflowing turbulent nature of the wind which, 
given the possible range of outflow, radial & turbulent ve- 
locities, may show absorption variability on timescales from 
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hours upwards. The two observation of Mrk 6 discussed here 
are each insufficiently long to probe absorption variation on 
timescales longer than ~4 hours making future long, or mon- 
itoring, observations of this source important to place any 
constraints on the nature of the wind. 



5 CONCLUSIONS 

A key feature of AGN unification schemes is the capability of 
absorbing material to change the appearance of a the central 
engine. To this end is it essential to understand the absorp- 
tion present in these systems both to confirm or refute the 
standard model and to allow us to make conclusions regard- 
ing the central engine and the primary emission regions in 
AGN. 

In this paper we have presented an analysis of the most 
recent XMM-Newton observation of Mrk 6. In particular, 
we test the hypothesis that the absorption can be well char- 
acterised by a warm absorption model that offers a more 
physical and testable description of the X-ray absorption 
than the standard, rather ad-hoc, partial covering model. 
We find that the simple warm absorption model, with only 
a single line of sight, provides an equally good statistical 
representation of the XMM-Newton CCD data as a partial 
covering model. Furthermore, we find that once the RGS 
data are included in the spectral fitting, the simple warm 
absorber model provided a very good fit to the data, with- 
out adding any additional components to the model, con- 
trasting with the partial covering model which requires the 
addition of either a low metalicity thermal plasma or black- 
body emission in order to provide a similar quality fit. 

Further support for the presence of warm absorption 
comes when the previous XMM-Newton observation is con- 
sidered. The absorbing column density has varied consider- 
ably between the two observations. Fitting the partial cov- 
ering model model to both sets of spectra does provide a 
reasonable fit and can explain the ~Nh changes but, crucially, 
the covering fractions of the two dominant line of sight in the 
model have remained the same. To put this in context, this 
implies a situation in which absorption in moving clumpy 
material relatively close to the source results in two distinct 
lines of sight to the source, both with separate absorption 
that varies considerably without resulting in any apprecia- 
ble change in the covering fractions involved. By contrast the 
warm absorber model needs only a second, higher column 
density but higher ionisation absorber to be present along 
the line of sight during the previous XMM-Newton obser- 
vation to explain the absorption variability. The presence 
of a second 'phase' of the warm absorption is simple to in- 
terpret; an origin for the absorbing material in an ionised, 
outflowing, clumpy accretion disk wind or skin of the molec- 
ular torus provide a natural explanation for the absorption 
variability. The densities implied by the ionisation state of 
both the original warm absorber and the second, higher col- 
umn density, higher ionisation, absorber are in good agree- 
ment with the densities found from simulations of these re- 
gions. Furthermore, the presence of several distinct ionisa- 
tion 'phases' is in good agreement with the results from the 
high resolution grating X-ray spectra of most bright type-1 
AGN. Observations of these sources often display complex 
absorption with several distinct, clearly identifiable, ioniza- 



tion phases (e.g. NGC 3783 - Netzer et al. 2003, Behar et 
al. 2003, NGC 3516 - Turner et al. 2005, etc.). 

Based on this evidence, coupled with the relatively ar- 
bitrary nature of the partial covering model, we strongly 
favour the interpretation of this absorption as originating in 
ionised clumpy material, probably resulting from an accre- 
tion disk wind. 
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